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'nI" I Standard QCD sum-rule analyses of the nucleon mass give results that are 

inconsistent with chiral perturbation theory due to an overly simple contin- 

■^ ' uum ansatz on the phenomenological side of the sum rule. We show that a 

careful treatment of the continuum, including vr-A^ states and other states 

with virtual pions, resolves the inconsistency associated with chiral logs. 
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While we believe that QCD is the theory underlying strong interactions, the problem 
of describing low-energy hadronic physics remains an active and controversial field. The 
essential difficulty is that QCD remains intractable in this regime; thus progress has been 
made primarily through partial treatments of the problem, where knowledge of QCD is 
supplemented by phenomenological input. Chiral perturbation theory |Q^ and QCD sum 
rules 0,^ are two largely orthogonal approaches of this variety that have played a significant 
role in describing and explaining low-energy hadronic phenomena. In this letter, we explore 
the relationship between the two approaches. 

It is generally believed that chiral perturbation theory accurately describes all low-energy 
observables of QCD in the limit of light current quark masses ||1],0]. In the present context, 
however, we will only be interested in the leading nonanalytic behavior of observables as 
a function of niq, the average of the up and down current quark masses. It should be 
noted that this nonanalytic behavior can be determined without the full machinery of chiral 
perturbation theory — it depends only on the existence of dispersion relations and a pseudo- 
Goldstone pion. 

The basic idea of QCD sum rules |^,Q is to extrapolate from the large spacelike momen- 
tum region, where we know how to treat QCD, down to the low timelike region relevant for 
hadronic physics. This requires three steps. First, a time-ordered correlation function of 
interpolating fields is calculated using the operator product expansion (OPE), which gives a 
large-momentum expression for the correlator; information about the nonperturbative QCD 
vacuum enters this description via nonvanishing condensates. Second, the correlator is re- 
lated to the spectral density using a dispersion relation; a model with a small number of 
parameters is then made for the spectral density, which gives a phenomenological description 
of the correlator. Finally, the OPE and phenomenological sides are matched in a manner 
based on the analytic properties of the correlator and asymptotic freedom, and the spectral 
parameters, such as masses and couplings, are extracted. 

In principle, an exact evaluation of the correlator in terms of QCD degrees of freedom 
must reproduce all of the low-energy hadronic physics, including the physics associated with 
chiral perturbation theory. However, all practical implementations of QCD sum rules are 
not exact, and it is not guaranteed that the relations thus obtained are consistent with 
low-energy constraints. 

In a previous work P, two of us (DKG and TDC) pointed out inconsistencies in the 
leading nonanalytic behavior in rriq between the usual QCD sum-rule treatment of the nu- 
cleon mass (based on a simple continuum ansatz in the spectral density) and the chiral 
perturbation theory description. The simplest QCD sum-rule formula for the nucleon mass 
1^ directly relates the nucleon mass to the quark condensate. The inconsistency originates 
from the fact that the quark condensate has a leading nonanalytic term proportional to 
TTiq In niq [^||,|ll , whereas the nucleon mass is known to have a leading nonanalytic term 
proportional to m^/^, but no mqlnniq term 0. In Ref. 0, it was argued that the origin of 
this problem clearly lies in the overly simple model for the continuum usually used in QCD 
sum rules. The consequence of this was nontrivial — a rough estimate of the uncertainty due 
to the inconsistent nonanalytic behavior was of order 100 MeV for the nucleon mass, with a 
comparable uncertainty for the a term. 

In this letter we will show that the inconsistencies between chiral perturbation theory 
and QCD sum rules disappear if the physics of virtual pions is properly taken into account 



on the phenomenological side of the sum rule. This is quite reminiscent of the case of octet 
baryons at finite temperature, as described in Ref. [0. The consequence of this for QCD 
sum-rule descriptions of the nucleon will be discussed at the conclusion of this letter. 

QCD sum-rule analyses of the nucleon mass Q are based on the time-ordered correlation 
function IIn{(i) defined by 

n^(g) = i /"(i^xe^«-^(vac|Tr/^(x)r7^(0)|vac) , (1) 

where |vac) is the physical nonperturbative vacuum state, and rj^ is an interpolating field 
with the spin and isospin of a nucleon, but with indefinite parity. We use the usual loffe 
interpolating fields for the proton and neutron P,pi]|: 



Vp = (^abc{u^C-f^Ub)j5j''dc , T]n = -€abc{dlCj^db)j5l^Uc , (2) 

where Ua and da are up and down quark fields (a is a color index), T denotes a transpose 
in Dirac space, and C is the charge-conjugation matrix. These choices of the interpolating 
fields are not unique; one can also choose 

Vp = (^abc{ulCa^^Ub)'j5(^'"'dc , r/^ = -eabcid'^Ca^^db)'j5(^'"'Uc , (3) 

but r]p and 7]n are the standard — and most effective — interpolating fields for QCD sum-rule 
studies of the nucleon |^,0. Our conclusions will not, however, depend on this choice. 
Lorentz covariance and parity invariance imply that the Dirac structure of IIn{(i) is of 



the form 11,12 



U^{q)^n,{q')+n,iq')i. (4) 

Asymptotic freedom and analyticity imply that the Lorentz scalar function Ili{s) {i = {1, q}) 
satisfies the sum rule [0,|1^,0 

/ dsW,{s)AU,{s) = dsWi{s)AUf''''{s) , (5) 

Jo Jo 

where the discontinuity, 

AUi{s) = lim Ui{s + ie) - Ui{s - ie) , (6) 

is proportional to the spectral density, nP^^(s) denotes the scalar function evaluated using 
the OPE, and the weighting function Wi{s) is an arbitrary entire function. The weight- 
ing function is chosen so as to improve the convergence of the OPE while simultaneously 
strongly weighting the nucleon pole contribution to the sum rule relative to the continuum 
contribution. The usual QCD sum rules based on the Borel transform 0J^ can be obtained 
by choosing 1^1(5) = Wq{s) = e~^^^^ , where M is known as the Borel mass. 

The correlator contains contributions from the nucleon pole and from higher-mass (con- 
tinuum) states with the quantum numbers of a nucleon; thus the correlator can be parame- 



terized as 11,12 



^Niq) = -\%j^j^ + IlTiQ) , (7) 

where Xn specifies the strength of the couphng between the interpolating field and the phys- 
ical nucleon state [see Eq. (pOl)], and 11™'^* (g) denotes the contribution from the continuum. 
On the other hand, the OPE expresses the correlator in terms of vacuum condensates. 
With the interpolating fields in Eq. (j^), one obtains [jBlpTSf 

n?^^(^) = ^ln(-.)(gg)-^^(gg)((^G)^) + --- , (8) 

nr(^) = -^M-^)-Y^in(-.)((^Gr)-^(gg)^ + --- , (9) 

where (qq) ~ — (225MeV)^ [qq = ^{uu + dd)] and {(gG)"^) ~ 0.5 GeV"^ are the quark and 
gluon condensates [3,H. 

The usual strategy for extracting the nucleon mass in QCD sum rules is to separate the 
nucleon pole contributions to the sum rules from the OPE and continuum contributions. 
For simplicity, we use the same weighting function W{s) in both sum rules. The nucleon 
mass is then obtained by taking the ratio of the sum rule for Hi to the sum rule for 11^ [0: 
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/ dtW{t)\AU^^^{t)- An^i°°* (t) 
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dsW{s) 



AIl?^^(s)-AUTHs 



(10) 



The formula for the nucleon mass in Eq. ([T0|) has a potential problem when approaching 
the chiral limit [^. As we will show below, AII^^^^ contains an 0{mq\nmg) chiral log, 
whereas AII^^^ and the nucleon mass [Q] lack this nonanalytic term near the chiral limit. 
This result is somewhat more general than that of Ref. |^, where it was pointed out that the 
leading term in AII^^^, which is proportional to (qq), has a well-known 0,|,|ll 0{mglnmq) 
chiral log. 

The physics associated with the chiral logs can be thought of as arising from virtual 
pions. As long as one studies the correlation function for g^ sufficiently far from the nucleon 
pole (Ig^ — M^l ^ 2Mjvm7r), one can use soft-pion theorems and phase-space arguments to 
show 

nr (g^) = (1 - 26)nr (g^) , (11) 

nf-(g^) = [l + 0(m^)]nf-(g^), (12) 

o o 

where we have defined m^ = — 2?Tig(gg)//^. All other quantities of the form x denote the 
value of X in the rriq = limit. We have also defined 

e^^^mlln^, (13) 

where Mq is an arbitrary constant. Changes in Mq can always be absorbed into changes 
in an analytic term proportional to m^; thus the substitution e —^ e + O(m^) is implicit. 



Equivalently, for any weighting function W{s) that has substantial strength over a region in 
s that goes well above the vr-A^ threshold 

/•OO /"OO o 

/ dsW{s)AU^''^{s)= rfsVr(s)(l-2e)AnOP^(s) , (14) 

Jo Jo 

/•OO rco o 

J ds ly (s) AnOP^(s) = j ds W{s) [1 + 0(m2 )] AnOPE(s) . (15) 

Note the restriction in Eqs. (p!T]) and (|12D to q^ far from the nucleon pole translates to a 
restriction in Eqs. (pjlf ) and (|l^) to weighting functions that cover a large range in s. 

In contrast, the nucleon mass is known to have no 0(m^ In m^) contributions. Near the 
chiral limit {ttIt^ -C Ma — Mjv), the nucleon mass is given by 

Mj, = Mj, + Aml-^^ml + -.. , (16) 

327r/2 

where A is an unknown constant . The essence of the difficulty is the following: Consider 
Eq. (^) with different values of the quark mass. The right-hand side of Eq. (jT^) has 
contributions to M^v of order mg In rUq coming from the quark-mass dependence of the OPE. 
On the other hand, we know that M^- has no term that goes as nig Inmg. 

In principle, there is no difficulty in reconciling the QCD sum-rule expression in Eq. (p!OD 
with the results of chiral perturbation theory. The continuum contribution in the numera- 
tor, for example, can have a nontrivial chiral behavior that precisely cancels the mglnrrig 
behavior. The difficulty is that, in most practical QCD sum-rule calculations, the model of 
the spectral function for the continuum is very crude: 

An^°'^*(s) = e{s - so)An°p^ . (17) 

The issue is whether there is any natural way for the continuum threshold sq in this model 
to vary with rUg in such a manner as to cancel the spurious mglnmg behavior. We will 
demonstrate that this is impossible. As we will show here, however, if a better model for the 
continuum is used, one that includes vr-A^ (and higher-mass) continuum states in a fashion 
consistent with chiral symmetry, then one automatically cancels the problematic rUg In mg 
behavior. 

The simplest derivation of Eqs. (|TT|) and (|12D is based on a spectral representation of the 
correlator,^ 

n,{g)=rd.>^i^^^^, (18) 

J-00 uj — qo 
where the spectral density is defined as 

PNiq) = (27r)3^(n|r/^(0)|m)(m|r7jv(0)|n)(5^(g-p„+p„) . (19) 



^We neglect imaginary infinitesimals in the denominator, since we are interested in the correlator 
at spacelike momentum transfers. 



The states \n) denote a complete set of energy eigenstates, which are also eigenstates of 
momentum due to translation invariance. The double sum in Eq. (19) is a convenient way 



to represent the contributions from positive- and negative-energy states; however, one must 
have either \n) = |vac) or \m) = |vac) in each term. To obtain the leading nonanalytic be- 
havior of IlAr, we consider N and vr-A^ intermediate states. It is clear that other intermediate 
states will not contribute to the leading nonanalytic term. 

We first evaluate the nucleon- and antinucleon-pole contributions to the correlator, which 
correspond to taking \n) = |vac), \m) = \N{k)) and \m) = |vac), \n) = \N{k)), respectively, 
in Eq. (|T^). Given the definitions 

(vac|r/;v(0)|iV(fc)) = XnUnW , (vac|r7^(0)|iV(fc)) = XNVN{k) , (20) 

where u^ik) and vj\f{k) are nucleon and antinucleon spinors, one readily obtains the result 
announced in Eq. (|^), 

from Eqs. (|TH|) and {^). _ 

Next we consider the contributions to the correlator from vr-A^ and vr-A^ continuum states. 
The corresponding states in Eq. (|19|) are of the form \n) = |vac), \m) = \TT{k')N{k)) and 
\m) = |vac), \n) = \TT{k')N{k)). To evaluate the resulting matrix elements, we exploit the 
transformation properties of the nucleon interpolating field under SU(2) axial rotations: 



where Qs is the axial charge. These transformation properties do not depend on the choice 
of interpolating field; the same relations hold for r/N — *■ tj'j^- Given the transformation 



properties for the interpolating field, one can derive the soft-pion theorem ||16 



(vac|C^|7r'^iV(A:)) = -^{ye.c\[QlCr,]\N{k)) + O (^^^] 

= ^75r'^(vac|C^|Ar(fc)) + O f ^) , (23) 

where Ahad is a typical hadronic scale assumed to be ~ 1 GeV. For g^ far from the nucleon 
pole, the soft part of the vr-A^ continuum yields the following contribution to IlN{q)'- 

TIn"" ''"'''{q) = -ei,Il^N%H ■ (24) 

The physical origin of this contribution is clear: The vr-A^ states become nearly degenerate 
with the A^ states, and the infrared behavior gives rise to the chiral log behavior. 

The second source of chiral logs is the A^ factor in the nucleon-pole contribution to the 
correlator. Consider the matrix element (vac|?77v|A^)- By treating the finite up and down 
quark masses as a perturbation, one obtains 



(n'^n'MN) = -_(vac|r/^|iV) + O { -^ ) . (28) 



(vac|r7^(0)|iV) = (vac|r/^(0)|iV)^^=o + ^J d^x {vac\T5jC{x)r]N{0)\N)c , (25) 

where we define 5C = —rriqiuu + dd). The leading nonanalytic behavior of the time-ordered 
product is given by two-pion intermediate states; the resulting matrix elements are approx- 
imated by taking the soft-pion limit. This yields 

(vac|?7Af|A^) = (vac|r7iv|A^)m,=o 

+ i™,(vac|™ + 3d|rt'> (.'■An.W I ^^^;j_L__ , (26) 

where all fields are evaluated at the origin. From PCAC one obtains 

m,(vac|MM + ddlir^ir'') = mid"'' , (27) 

and the matrix elment of t]^ is obtained by applying the soft-pion theorem twice: 

— (vac|r/^|iV) + -^ 

One then obtains the result 

{Ya.c\r]N\N) = (l - (vac|r7^|Ar)^^=o , (29) 

o 

which implies Aat = (1 — e/2)XN. Furthermore, the nucleon-pole contribution to the corre- 
lator then behaves as 

nr(?) = (i-e)nr(?)- (30) 

We must also determine the chiral behavior of the continuum contributions to the cor- 
relator. The effect of vr-A^ continuum states on U^ has been discussed above [see Eq.(p^)]. 
The chiral behavior of the remainder of the continuum can be deduced as follows: In the 
chiral limit, this part of the continuum can be represented to arbitrary accuracy by a dis- 
crete set of poles, A^*, A^**, . . ., of positive or negative parity. Away from the chiral limit, 
each A^* is dressed with pions in the same way as the nucleon; thus the leading nonanalytic 
contributions to this part of the correlator arise from the chiral expansion of Aat* and from 
n-N* states. One then obtains relations analogous to Eqs. ( PPD and ( pij ) for both positive- 
and negative-parity states. The expansion of the continuum contribution to IItv about its 
value in the chiral limit, for g^ far from the nucleon pole, is thus given by 

nr*(?) = -e75n°''(g)75 + (1 - 6)nr (?) - n^iiTi^H ■ (3i) 

The first term in the preceding equation arises from the fact that, while the nucleon pole 
term is, by construction, not part of the continuum, the vr-A^ states that give the chiral log 
in Eq. p^ ) are part of the continuum. 

By combining Eqs. ( pOD and (|3TD , we can determine the chiral expansion of the correlator. 



Decomposing IIn{(i) into scalar functions according to Eq. (|), we obtain 



ni(g2) = (l_2e)ni(g2), (32) 

n,(g^) = [l + 0(rn^)]n,(g^). (33) 

These results can be used at high spacehke momentum transfer, where they must match the 
OPE description. Thus one obtains the resuhs stated in Eqs. ( pJ]) and (0): 

n?P^(g2) = (1 - 2e)n?P^(g2) , (34) 

nr"(g^) = [l + 0(^^)]nf^(g^). (35) 

Since Eq. ( |34D holds for all spacelike g^, it follows that all terms in Ilf^^ must contain 
the (1 — 2e) factor. In Ref. |^, it was demonstrated that the leading term in 11^^^, which 
is proportional to {qq), behaves this way; here we see that all of the other terms do too. 
Similarly, all terms in U'?^^ must go as [1 + Olm^)]. 

Let us now use the known chiral expansions of 11^^^ and 11™"* to see whether there is a 
chiral-log contribution to M^ in Eq. (|T0|). Inserting the results of Eqs. (|34D, (|35D , and (|3TD 
into Eq. (^) gives 

Mn = Mn + 0{ml) . (36) 

Thus we see that, if the continuum contribution to the correlator is treated carefully, the 
spurious chiral-log contributions to the nucleon mass, discussed in Ref. P, cancel. Thus 
the QCD sum-rule prediction is consistent with the chiral perturbation theory description 
of the nucleon mass given in Eq. (|T^ to 0{mq). However, it lacks the 0{my^) contribution, 
which must arise from subleading nonanalytic contributions to the nucleon correlator. 

While it is clear that the inclusion of virtual pions automatically removes the spurious 
chiral logs, it is worth exploring the question of whether the usual continuum ansatz in 
Eq. {\[7\] can accomplish the same results. This is impossible. Decomposing the continuum 
contribution to the correlator given in Eq. (^) according to Eq. (^), one obtains 

w^-^'iq^) = -enr''(g') + (1 - 2e)nr*(g') , (37) 

Ul°^\q') = elif^q') + [1 + 0{^lWr\q') ' (38) 



and it is straightforward to see that the form given in Eq. ([I7|) is inconsistent with these 
results. Physically, one can understand this by noting that there is spectral strength very 
near the nucleon pole coming from low-momentum vr-A^ states, which give rise to part of the 
chiral log. Alterations in the value of the continuum threshold, however, can only change 
the spectral density far above the nucleon pole. 

Finally, we consider the implications of our results for QCD sum rules for the nucleon. It 
seems clear that one should include low-lying tt-N continuum states on the phenomenological 
side of the sum rule; this will remove the spurious chiral logs and remove one source of 
uncertainty. The uncertainty from this effect has been estimated to be ~ 100 MeV for 
both the nucleon mass and a term 0. However, we note that, although this removes the 
rather unsatisfactory chiral log behavior, it does not completely remove the uncertainties. 
In particular, chiral perturbation theory gives a term proportional to m^^^ as the leading 



nonanalytic contribution to the nucleon mass, which arises from vr-A^ intermediate states in 
the nucleon propagator. We expect this term to originate in a similar manner in a QCD 
sum-rule analysis. It seems certain that the vr?!'^ term is lost in the present QCD sum-rule 
description due to the use of the soft-pion theorem in evaluating matrix elements involving 
T[-N intermediate states [see Eq. (|23|)]; presumably, it can be reproduced with an adequate 
description of the soft-to-hard pion correction. In fact, the w?!'^ term is recovered using a 
simple model for the soft-to-hard correction based on pseudovector vr-A^ coupling; however, 
we know of no satisfactory, model-independent approach. The lack of the vr?!"^ term leads 
to uncertainties in the sum-rule results of ~ 15 MeV and ~ 20 MeV for the nucleon mass 
and 0" term, respectively 0. 
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